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Abstract 


We perform a Hamiltonian reduction of spherically symmetric Einstein grav¬ 
ity with a thin dust shell of positive rest mass. Three spatial topologies are 
considered: Euclidean (M^), Kruskal (F^ x M), and the spatial topology of a 
diametrically identihed Kruskal (MP^\{a point at inhnity}). Eor the Kruskal 
and MP^ topologies the reduced phase space is four-dimensional, with one 
canonical pair associated with the shell and the other with the geometry; 
the latter pair disappears if one prescribes the value of the Schwarzschild 
mass at an asymptopia or at a throat. Eor the Euclidean topology the re¬ 
duced phase space is necessarily two-dimensional, with only the canonical 
pair associated with the shell surviving. A time-reparametrization on a two- 
dimensional phase space is introduced and used to bring the shell Hamilto- 
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nians to a simpler (and known) form associated with the proper time of the 
shell. An alternative reparametrization yields a square-root Hamiltonian that 
generalizes the Hamiltonian of a test shell in Minkowski space with respect 
to Minkowski time. Quantization is briefly discussed. The discrete mass 
spectrum that characterizes natural minisuperspace quantizations of vacuum 
wormholes and MP^-geons appears to persist as the geometrical part of the 
mass spectrum when the additional matter degree of freedom is added. 

Pacs: 04.20.Fy, 04.40.Nr, 04.60.Kz, 04.70.Dy 
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I. INTRODUCTION 


In classical general relativity, every three-manifold occurs as the spatial topology of 
a globally hyperbolic vacuum spacetime. In a canonical approach to quantum gravity, the 
spatial topology is frozen, and one can ask for ground states corresponding to each topology.[| 
Spherically symmetric minisuperspaces provide simple models for the quantization of 
geometries with non-Euclidean topology. The spatial topologies consistent with spherical 
symmetry and asymptotic flatness are the wormhole x M of the Kruskal geometry 
with two asymptopias, and the geon, a manifold with a single asymptopia obtained by 
removing a point from the compact manifold MP^. This last manifold is the space acquired 
from Kruskal geometry by identifying diametrically opposite points on an U -|- U = constant 
slice, with U and V the usual Kruskal null coordinates |^. 

A reduced phase space formalism for spherically symmetric vacuum Einstein gravity in 
four spacetime dimensions has been considered by several authors &II3-0 In the present 
paper we add to spherically symmetric Einstein gravity an idealized, inhnitesimally thin dust 
shell of positive rest mass. The equations of motion for such a shell follow easily from Israel’s 
junction condition formalism p2|-|35[| , and a number of workers have proposed actions from 
which these equations can be derived • Our main purpose is to find an action for this 


system by an explicit Hamiltonian reduction, treating both the geometry and the shell as 
dynamical, and retaining the full dynamics allowed by the choice of the spatial topology. 
Two issues require particular care. First, as general relativity is a nonlinear theory, 

H 


introducing a distributional source faces well known subtleties 


The special case of a 


source concentrated on a hypersurface of codimension 1 is fortunate, as Einstein’s equations 
can then be given an unambiguous distributional interpretation, and this interpretation 
reduces to Israel’s junction conditions when the source is a pure delta-function on the surface 
0- However, we wish to go further and write an action principle from which the held 
equations would arise as variational equations. In such an action principle one presumably 
needs to be able to vary the action with respect to both the metric and the shell variables, 
with the variations remaining independent in some suitable sense. This brings in not only 
the regularity properties of the spacetime at the shell, but also the regularity properties of 
the coordinates in which the action is written. 

We will not hnd an action principle whose variational equations would be fully distribu- 


^Even in a theory that permits topology change, topologies threaded by electric or magnetic flux 
in source-free Einstein-Maxwell theory (or in higher-dimensional gravity with Kaluza-Klein asymp¬ 
totic behavior) cannot evolve to Euclidean space, if they have a net asymptotic charge. If there is 
a nonsingular quantum theory of such a system, it must allow a ground state with nonzero asymp¬ 
totic charge and non-Euclidean topology. Topological geons with half-integral angular momentum 
in a quantum theory of gravity would similarly be unable to settle down to Euclidean topology. 


^Eor extensions to related theories, including spherically symmetric Einstein-Maxwell theory and 
lower-dimensional dilatonic theories, see Refs. ||l^,13-^]. Eor discussions within the Euclidean 
context, see for example Refs. |22-31] and the references therein. 
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tionally consistent at the shell. However, the ambiguity in our variational equations will be 
localized into the single equation that results from varying the action with respect to the 
shell position coordinate. When the ambiguous contribution to this equation is interpreted 
as the average of its values on the two sides of the shell, as is necessitated by consistency with 
the rest of the equations, we correctly reproduce the content of Israel’s junction condition 
formalism. At a somewhat formal level, our action will be manifestly invariant under the 
Hamiltonian version of spacetime coordinate transformations preserving spherical symmetry. 

Second, one needs to choose the falloff and boundary conditions at the asymptopias. We 
shall set the asymptotic momenta to zero, but the values of the Schwarzschild masses at 
the asymptopias will be left free to emerge from the dynamics. Our spacelike hypersurfaces 
will not be asymptotic to hypersurfaces of constant Minkowski time, but the foliation is 
nevertheless asymptotically Minkowski in the relevant sense. In particular, the generator of 
unit time translations at the infinity is the Schwarzschild mass. 

We shall find that the reduced phase space is four-dimensional with the Kruskal and MP^ 
topologies, and two-dimensional with the topology. With each topology, one canonical 
pair is associated with the shell motion, but with the Kruskal and MP^ topologies there is 
also a second canonical pair, associated with the dynamics of the geometry. In the limit 
where the shell is removed, this reproduces results previously obtained in the Hamiltonian 
vacuum theories with the Kruskal and RP^ topologies [^,0. 

For the non-Euclidean topologies, the canonical pair associated with the geometry dis¬ 
appears if one prescribes by hand the mass at one infinity in the Kruskal topology, and 
the mass at the wormhole throat in the RP^ topology. All three reduced phase spaces then 
become two-dimensional, and they can be treated on an essentially equal footing. 

We next introduce a formalism for reparametrizing time in a Hamiltonian theory with a 
two-dimensional phase space. Applying this formalism to our two-dimensional phase spaces, 
we redehne the coordinate time to coincide with the proper time of the shell and thereby 
obtain a Hamiltonian that can be given in terms of elementary functions. This Hamiltonian 
is known [^, but the fact that it emerges from a minisuperspace framework is new. An 


alternative choice for the coordinate time yields a Hamiltonian that generalizes to our self- 
gravitating shell the familiar Hamiltonian of a spherical test shell in Minkowski 

space. 

The paper concludes with a discussion of the prospects for quantization. Quantization 
of the vacuum case is revisited to emphasize choices that lead to discrete or continuous mass 
spectra. The additional degree of freedom provided by the shell does not appear to qualita¬ 
tively alter these choices. Like the Jain-Schechter-Sorkin quantum-stabilized Skyrmion |[43|| , 


the minisuperspace geons provide an example of field configurations that have quantum, but 
not classical ground states; both are held theory analogues of the quantum-stabilization of 
the hydrogen atom. Whether the ground state of geons is an artifact of the reduction of the 
degrees of freedom is, of course, an open question, but the geometrical ground state appears 
to persist when the shell’s degree of freedom is added. 

With one asymptotic or interior mass hxed, the implicit Hamiltonian we obtain prior to 
time-reparametrization was found by Kraus and Wilczek in the limit of a massless 

shell, and it could easily be found from what they present also for the massive case. A related 
reduction technique was used earlier by Fischler et al. in a minisuperspace treatment 


of a bubble wall, and recently generalized by Kolitch and Eardley [^|. For a hat geometry 
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interior to the shell, our proper-time Hamiltonian has been considered classically in Ref. 
and quantum mechanically in Ref. |^]. For a flat geometry interior to the shell, quantization 


using the square root Hamiltonian has been considered in Refs. 


Latin tensor indices a,b,... indicate abstract spacetime indices, 
units, h = c = G = 1. 


We work in Planck 


II. HAMILTONIAN FORMULATION FOR SPHERICALLY SYMMETRIC 
GEOMETRY WITH A DUST SHELL 

In this section we set up a Hamiltonian formulation for spherically symmetric gravity 
coupled to a thin dust shell. We pay special attention to the smoothness of the gravitational 
variables and to the global boundary conditions. 


A. Bulk action 

A spherically symmetric spacetime metric can be locally written in the Arnowitt-Deser- 
Misner (ADM) form 


ds'^ = -N^dt^ + A\dr + N^dt)^ + R^dn‘^ , ( 2 . 1 ) 

where is the metric on the unit two-sphere, and N, N'', A, and R are functions of t 
and r. Issues of smoothness and global structure will be addressed below. We denote the 
derivative with respect to t by overdot, and the derivative with respect to r by prime. 

The matter consists of a thin shell of dust, with a hxed positive rest mass m. We write 
the trajectory of the shell as r = t{t). Denoting by N{t), N^{t), A{t), and R{t) the values 
of N, A, and i? at r = r, 

R{t) := R(t,x{t)) etc, (2.2) 


the Hamiltonian action for the shell is 


5'sheii = / dt [pi- NR p2A-2 + m'^ + N''p 


(2.3) 


with p being the momentum conjugate to r. One can think of the shell as a spherically 
symmetric cloud of massive relativistic point particles. 

The Lagrangian gravitational action for the geometry m is obtained by integrating 
the Lagrangian density (^R — Kab + over the two-sphere p|,0,||, 

After A and R are replaced by their conjugate momenta. 


44,46,48 


^r = ~(R-N'R') 


R 

N 


A - (N^A)' 


(2.4a) 

(2.4b) 


the Hamiltonian bulk action for the coupled system reads 
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(2.5) 


= dt 


pr + dr { ttaA + nnR - NH - 


where the super-Hamiltonian Ti and the supermomentum are given by 

_ Avr^ ttattr RR" RR'R' R'^ a 

+ \J p2A“2 frii — x) , 

R-r = t^rR! — ttaA — p(5(r — r) . 


(2.6a) 

(2.6b) 


We shall hrst discuss the smoothness of the gravitational variables, and then the boundary 
terms to be added to the bulk action. 


B. Smoothness 

In the presence of a smooth matter distribution, one can assume the spacetime metric 
to be smooth ((7°°). In the idealized case of an inhnitesimally thin shell, the metric can be 
chosen to be continuous but not, in general, differentiable across the shell p2|-[35[|. In the 
particular case of a spherically symmetric dust shell, Einstein’s equations imply that the 
extrinsic curvature of the shell history is discontinuous both in its angular components and 
in its component along the shell four-velocity. If the metric is taken continuous, we must 
therefore accommodate discontinuities in R' and in at least someQ of A', N\ and (A^'")'. We 
would like both the action ( p.5|) and its local variations to be well dehned, and such that 
the resulting variational equations are equivalent to Einstein’s equations with a dust shell. 

To proceed, we assume that the gravitational variables are smooth functions of r, with 
the exception that N', A', R', tta, and ttr may have hnite discontinuities at isolated 

values of r, and that the coordinate loci of the discontinuities may be smooth functions of t. 
It will be shown that the resulting variational principle is satisfactory in the above sense, 
provided one of the variational equations is interpreted as the average of a discontinuous 
quantity over the two sides of the shell.|^ All the terms under the r-integral in the action 
( p.5|) are well dehned in the distributional sense. The most singular contributions are the 
explicit matter delta-contributions in the constraints, and the implicit delta-functions in R" 
and tta. All these delta-functions are multiplied by continuous functions of r. The remaining 
terms are at worst discontinuous in r. The action is therefore well dehned. 


^By continuity of the metric, R{t) is well defined for all t. Taking the total time derivative of 
(|2.2| ) shows that AR = —xAR', where A denotes the discontinuity across the shell. Similarly for 
A, N, and N'". Continuity of A', N', and (N'")' would therefore imply that the extrinsic curvature 
of the shell history is discontinuous only in its angular components. 

'^Because the constraint equations enforce smoothness of the metric outside the shell, our differ¬ 
entiability assumptions can probably be relaxed. 
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Local independent variations of the action with respect to the gravitational and matter 
variables give the constraint equations 


Lf = 0 , 
Lf, = 0 , 


and the dynamical equations 


R = 


TTa = 


r = 


P = 


iV( 

V R2 r) 

+N'^R! , 
R 

N 

vri (R' 

y 

R2 ^A 

N 

Avr^ 

R3 R2 


Np 

}R‘\J p2A“2 -p rn? 


7VA'p2 


y) + , 


2p^h(r — r) 


+ 1 + 

K.^\J p2A-2 -)- 777,2 

R'\'^ fN'R\' 

A 


N'RR' 

A2 




A 


) +(A'’'irR)' . 




N'Rp‘^K-'^ + w? + p{N^y . 


(2.7a) 

(2.7b) 


(2.8a) 

(2.8b) 

(2.8c) 

(2.8d) 

(2.8e) 

(2.8f) 


A^y p2A“2 + rrR 


With the exception of (|2.8f|) , all the equations (p.7|) and (|2.8|) are well dehned in a distribu¬ 
tional sense.[| What needs to be examined is the consistency and dynamical content of the 
well-dehned equations, and the interpretation of the single troublesome equation (|2.81|) . 

As a preliminary, consider the variation of the matter action S'sheii ( |2.3| ) with respect to 
the metric. From the definition of the stress-energy tensor. 






(2.9) 


and the equation of motion (|2.8e|) , we find that the surface stress-energy tensor of the shell 
[Ref. 1^, equation (21.163)] takes the form 


gab ^ 


m 


CL^.h 


47rR2 


■ u u 


( 2 . 10 ) 


®The constraint equations (^) contain explicit delta-functions in r from the matter contribution 
and implicit delta-functions in R" and vr)^. The right-hand sides of (|2.8a|) and (2.8b) contain at worst 


finite discontinuities, and the right-hand sides of ( 2.8c ) and ( 2.8d[) contain at worst delta-functions. 
This is consistent with the left-hand sides of ( |2.8a| )-( ^.8d ), because the loci of nonsmoothness in 
A, R, TTA and tttj may evolve smoothly in t. Note that both the explicit matter delta-functions and 
the implicit delta-functions in R" and are multiplied by continuous functions of r. 
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where is the four-velocity of the shell, normalized in the usual way u'^Ua = —1. This 
conhrms that the shell indeed consists of pressureless dust, with surface energy density 
m/(47ri?^) and total rest mass m. 

Also, recall that the full content of the Einstein equations at the shell is encoded in Israel’s 

We shall refer to the two sides of the shell as the “right-hand 
left-hand side”, in view of the Penrose diagram in which two partial Kruskal 


junction conditions p^ , |33 
side' 


and the 

diagrams are joined to each other along the shell trajectory. Israel’s junction conditions then 
read 


- 87i{Sab - iKbS) = K+ - K 


^ab 


( 2 . 11 ) 


where Ua is the right-pointing unit normal to the shell history, hab = Qab — naUb is the 
projector to this history, Kab = h'^ah’^b^end is the extrinsic curvature tensor, and the signs 
± refer respectively to the right and left sides of the shell. With ( p.lOl) , and with Kruskal 
geometries of masses M± on the two sides of the shell, the angular components of (p.ll|) 
read 


m 


\ 



+ 1 - 


2M, 


R 


— e_ 





+ 1 - 


2M_ 


R 


( 2 . 12 ) 


where r is the shell’s proper time. e+ = 1 (e_ = 1) if, when viewed from the geometry 
right (left) of the shell, the shell is in the right-hand-side exterior region of the Kruskal 
diagram, or if the shell is in the white-hole region and moving to the right, or if the shell is 
in the black-hole region and moving to the left. Otherwise e+ = — 1 (e_ = — 1). It can be 
verihed that the shell motion is completely determined by the single equation ( p.l2|) and the 
vacuum Einstein equations away from the shell. In particular, these equations imply that 
the tangential component of ( |2.11| ) is satished. A more explicit discussion can be found in 
Ref. m. 


Now, away from the shell, equations (E3) and (|2.^) are well known to be equivalent to 


Einstein’s equations. At the shell, the constraints (p.7|) read 


AR' = 




R 


AtTa = -y , 
A 


where we have adopted the notation (cf. footnote 1) 

A/;= hm [/(r + e)-/(t-e)] 

e^0+ 


(2.13a) 

(2.13b) 

(2.14) 


identifying r > r (r < r) as the right (left) side of the shell. Using (|2.8e| ) and ( |2.1CI| ), one 
hnds that (|2.13a|) is equivalent to (p.l2|) . Our Hamiltonian equations away from the shell 
and the Hamiltonian constraint ( p.7a| ) at the shell therefore form a system that is equivalent 
to the correct dynamics for the shell. When these equations hold, it can be verihed that 
equation (|2.13bD is proportional to (|2.13a|) by the nonsingular factor .R(r-|- N^)/N, and the 


































momentum constraint thus contains no new information. Similarly, it can be verified that 
the delta-parts in (p.8c|) and (p.8d|) reduce to identities and contain no new information. 
Finally, (|2.8a| ) and (p.8b| ) contain no delta-functions, and thus no new information, at the 

shell. 

The single remaining equation of motion is ( p.81|) . If the geometry were smooth at 
the shell, equations (p.8e|) and (|2.8fl) would by construction be equivalent to the geodesic 
equation for the shell, as can indeed be explicitly verihed. If the ambiguous spatial derivative 
terms in ( p.81]) are evaluated on the left (right) side of the shell, (|2.81]) thus implies the 
geodesic equation for the shell in the geometry on the left (right). However, these two 
geodesic equations are mutually inconsistent, and the shell motion implied by the rest of the 
equations is not geodesic in either of the two geometries. Instead, the rest of the equations 
imply that the left-hand side of (|2.81|) is equal to the average of the right-hand side over the 


two sides of the shell. (For the generalization of this observation to nonspherical dust shells, 
see Exercise 21.26 in Ref. ||^.) Therefore, if the ill-dehned right-hand side of ( |2.81|) is given 
this average interpretation, our Hamiltonian formalism reproduces Einstein’s equations for 
the dust shell. 

We are not aware of an a priori justihcation of the averaged interpretation of the right- 
hand side of ( ^.81|) . This interpretation is merely forced on us by the rest of the variational 
equations. In a strict sense, we therefore regard the variational principle as inconsistent, 
and the averaged interpretation of (|2.8f|) as put in by hand. Nevertheless, we shall proceed 


with this variational principle. It will be seen in section |IT| that the Hamiltonian reduction 
can be carried through with no apparent inconsistency. 

One check on the consistency of the formalism is that the Poisson brackets of our con¬ 
straints can be verihed to obey the radial hypersurface deformation algebra [^], as in the 
absence of the shell. If we denote by A/'(r) and M^{r) smooth smearing functions of compact 
support, the algebra has the form 


dr 

j drN'2'H 

^ drU^Hr, 

J drMH 

drAflHr, j 

drMyUr 


= / , 


= / dr -M',{M\y]nr 


(2.15a) 

(2.15b) 

(2.15c) 


C. Asymptopias and boundary terms 

We now turn to the global properties of the geometry. In this section we take the spatial 
topology to be that of the extended Schwarzschild geometry, x M = S'^\{two points}, the 
omitted points being associated with asymptotically hat asymptopias. The spatial topologies 
MP^\{a point at inhnity} and will be discussed respectively in sections |V]and lYJ 

At a general level, restricting the asymptotic behavior of an asymptotically hat system 
allows one to hx the momentum, angular momentum, and mass at spatial inhnity. In a 
quantum theoretic context, to restrict in this way the asymptotic behavior of the operator 
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and its conjugate momentum 7r“^ is equivalent to restricting the state space to an eigen- 
subspace of fixed total momentum, angular momentum, or mass. In our particular case of 
spherical symmetry, the angular momentum is necessarily zero. It would be consistent with 
spherical symmetry to allow a nonzero momentum at inhnity (in the classical framework, 
this would mean allowing boosted Schwarzschild solutions), but for our purposes this free¬ 
dom does not appear signihcant, and we shall set the momentum at inhnity to zero. We 


shall, however, retain the freedom associated with the system’s total mass. 

We take the coordinate r to have the range —cxo < r < oo. At the asymptopias r —> ±cxo, 
we introduce the falloff 

A(t,r) = l + , (2.16a) 

R{t,r) = \r\+0°°(^\r\~^-^^ , (2.16b) 

nA{t,r) = ^/2M±\r\ + , (2.16c) 

7rR(t,r) = y^ + 0°°(|rr^“^) , (2.16d) 

Ar(t,r) = l + 0°°(^|r|"^) , (2.16e) 

Ar"(t,r) = + , (2.16f) 


where M±{t) are positive-valued functions of t, and /3 is a positive parameter that can be 
chosen at will. indicates a quantity that is bounded at inhnity by a constant times its 
argument, with the corresponding behavior for its derivatives. 

It is straightforward to verify that the falloh (p.l6|) is consistent with the constraints and 
preserved in time by the dynamical equations. When the equations of motion hold, are 
independent of t, and their values are just the Schwarzschild masses at the two asymptopias. 
Using (p.l2 ), it is easy to show that the existence of two asymptotically hat inhnities implies 
that both asymptotic Schwarzschild masses in the classical solutions are necessarily positive. 
The assumption M±{t) > 0 in (|2.16|) does therefore not exclude any solutions. 

The falloh ( p.l6|) is not consistent with the conventional fallohs (see, for example. Refs. 
i6|j50[|) in which the hypersurfaces of constant t are asymptotic to hypersurfaces of constant 
Killing time when the equations of motion hold. Instead, the falloh ( |2.16|) is asymptotic to 
the ingoing spatially hat coordinates |5T|-^ , individually near each asymptopia. When M± 
are constants and all the 0°°-terms vanish, (^.16 ) yields the Schwarzschild metric in the 
ingoing spatially hat coordinates, separately for r > 0 and r < 0. Our reason for adopting 
( p.l6| ) is that the spatially hat gauge will prove useful in the Hamiltonian reduction in section 

EH [0- 

In a variational principle that does not hx the values of M±, the bulk action ( p.5|) must 
be amended by a boundary action. With our falloh (|2.16|) , the spatial metric approaches 
hatness at r —> ±oo so fast that the variations of R and A give rise to no boundary terms 
from the inhnities. The only nontrivial boundary term arises from integrating by parts the 
term J dt J dr K{5 tta)'^ associated with the momentum constraint. This boundary term 
is canceled if we add to the bulk action (p^.5|) the boundary action 
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SdJ: — — 


dt (M+ + M_) 


(2.17) 


The generator of unit time translations at the inhnities is therefore still the Schwarzschild 
mass, despite the unconventional falloff. 


III. REDUCED PHASE SPACE FORMULATION 

In the absence of the shell, the Hamiltonian reduction of our theory with a technically 
different but qualitatively similar falloff at the two asymptopias was discussed in Ref. p. 
When the asymptotic masses are not fixed, it was found that the reduced phase space is two- 
dimensional, whereas if one asymptotic mass is fixed, the reduced phase space has dimension 
zero. As the shell brings in one new canonical pair but no new constraints, one expects that 
the reduced phase space of our theory is four-dimensional when the asymptotic masses are 
not fixed, and two-dimensional if one asymptotic mass is fixed. In this section we shall verify 
this expectation by an explicit Hamiltonian reduction. 


A. Gauge transformations and the Hamiltonian reduction formalism 


In the Hamiltonian theory formulated in section ||, the variables (A, R, tta, tt/j, r, p) con¬ 
stitute a canonical chart on the phase space S, while N and act as Lagrange multipliers 
enforcing the constraints. As the Poisson bracket algebra (p.l5|) of the constraints closes, 
we have a hrst class constrained system . 

Let r denote the constraint hypersurface (|2.7|) in S. We take gauge transformations 
to mean the transformations on P generated by the constraints.^ Denoting the smearing 
functions by Af{r) and as in (|2.15|) , the smeared Hamiltonian constraint transforms 

an initial data set (A, R, tta, r, p) G P by the time evolution associated with TV, and the 
smeared momentum constraint transforms the initial data set by the spatial diffeomorphism 
associated with M''. The smearing functions must fall off so fast that the transformations 
become trivial at the infinities and the falloff ( p.l6|) is preserved.]] 

By definition, the reduced phase space f consists of the equivalence classes in P under 
gauge transformations. The symplectic form uj on 5, 


a; := 5p A (5r -|- / dr {6 tir A <5A -|- 6 tir A 6 R) 


(3,1) 


induces a symplectic form a> on P. Here, and from now on, <5 denotes the exterior derivative 
on the (functional) spaces in question. 


^See, for example. Ref. |^|. Note that this is distinct from, although closely related to, the gauge 


transformations that act on the histories on which the action is defined |M,56 


^One could consider an extended phase space that contains N and A'’ as new coordinates and 
their conjugates ttat and iiRir as new momenta. We shall, however, not need this extension. 
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We wish to implement this Hamiltonian reduction, hnding a> in an explicit symplectic 
chart on f. Our implementation will consist of the following three steps: 

1. Consider hrst F. At the shell, we have already seen that the full content of the 
constraints is encoded in equations (|2.13|) . Away from the shell, the constraints can be 
solved explicitly for the gravitational momenta as 


tta = R^{R'/kf - 1 + 2M±/R , 

_ A [{R/k){R'/ky + {R'/kf - 1 + M±/R] 
~ - 1+2M±/i? 


(3.2a) 

(3.2b) 


with the upper (lower) signs holding respectively for r > r (r < r). We have chosen the sign 
of the square root in ( ^.2|) so as to agree with the falloff (|2.1ti|) . This choice will lead to a 
reduction that will cover the black hole interior but not the white hole interior. 

2. To pass from T to F, we choose a gauge: we specify in F a hypersurface H that 
is transversal to the gauge orbits, so that each point in (an open subset of) f has a 
unique representative in H. This dehnes an isomorphism between H and (the open subset 
of) F. In order to choose the gauge in practice, we note that away from the shell, a point 
(A, i?, tta, tt/j, r, p) G F is an initial data set for the vacuum Einstein equations with spherical 
symmetry. Any vacuum initial data set has a unique time evolution, and, by Birkhoff’s the¬ 
orem, the resulting subspacetimes left and right of the shell are isometric to regions of two 
Kruskal spacetimes with the respective masses M_ and M_|_. A solution to the constraint 
equations can thus be regarded as two parametrized partial spacelike hypersurfaces in the 
two Kruskal spacetimes, joining appropriately at the shell. In this picture, a gauge choice 
means making a particular choice for these two partial spacelike hypersurfaces in the two 
Kruskal spacetimes, in a way that joins appropriately at the shell and is compatible with 
the falloff at the inhnities. 

3. To hnd the symplectic form cD on F, it is convenient hrst to hnd the corresponding 
Liouville form. Recall that on S, the Liouville form corresponding to our canonical chart 
(A, R,7rA,7r/{,r,p) is 


6 := p(5r -I- 


dr (tta^A -|- ttr6 R) 


(3.3) 


Pulling 6 back to H yields on H the Liouville form 9r, and cDjj := 69r is the symplectic 
form on on H that corresponds to a) on (the isomorphic open subset of) F. 

In view of the description of the gauge choice in step 2, a technical point in step 3 arises 
from the fact that although M± are constants in the time evolution of a given initial data 
set, they are not constants as functions on Ff, and their exterior derivatives may contribute 
to the pullback of 9 (|3.3|) . Put differently, a generic path in f need not correspond to a 
partial foliation of a single Kruskal geometry on either side of the shell. 

To complete steps 2 and 3, we need to specify the gauge. This will be described next. 


B. Gauge choice 

Our gauge choice involves taking the intrinsic metric on the spacelike hypersurface to 
be flat, with the exception of certain transition regions that are eventually taken to be 
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vanishingly narrow. The possible locations for the transition regions depend on whether the 
shell trajectory is visible to the right-hand-side future null infinity, the left-hand-side future 
null infinity, or neither.^ We now make the simplifying assumption that part of the shell 
trajectory is visible to one future null infinity, and we take this infinity to be on the right. 
This is arguably the situation of physical interest for an observer in the asymptotically flat 
region. 

Thus, fix an initial data set (a point in T), and consider the classical spacetime that 
is its time evolution. We assume that in this spacetime, the shell trajectory intersects the 
right-hand-side exterior region in the Kruskal geometry right of the shell. The shell equation 


of motion ( | 2 . 12 | ) then implies M_|_ > M_, and the trajectory intersects the right-hand-side 
exterior region also in the Kruskal geometry left of the shell. It follows that e_ = 1 on all of 
the trajectory, whereas e+ = 1 when R is sufficiently large (in particular, when R > 2M_|_) 
but e+ = —1 as i? — 0 . 

On this spacetime, we introduce two local charts, Ci and O 2 , as follows: 

Suppressing the angles, let the coordinates in the chart Ci be (ti,ri), with ri > 0. The 
metric reads 


= —dt\ -I- dri + 


ds^ = —dt\ -|- I dri + 


2M_ 

n 

2 M 4 

n 


dti + r\dVP‘, 0 < ri < r — / 


dti -I- rldQ"^, X < ri 


(3.4a) 


(3.4b) 


where I is a positive parameter. The two metrics shown in ( p.4|) are the ingoing right- 
hand-side spatially flat charts in Kruskal manifolds with the respective masses M_ and M+ 
If taken individually for 0 < ri < cx) and —00 <ti<oo, each of these two metrics 


would cover the upper right half (that is, the right-hand-side exterior and the black hole 
interior) in the respective full Kruskal diagrams. With the domains indicated in (p.4|) , the 
combined chart is spatially flat with mass M_ for ri < x — I, and spatially flat with mass 
M 4 for ri > X. The chart in the transition region x — I < ri < x will be specified below. 

Let the coordinates in the chart C 2 be {t 2 ,r 2 ), with r 2 < 0. The metric reads 


= —dti 


+ -dr2 + 


'2M_ 

1^2! 


dig 


-I- rldVt^, r 2 < 0 


(3.5) 


We identify C 2 as the ingoing left-hand-side spatially flat chart in a Kruskal manifold with 
mass M_, with r 2 —>■ —00 giving the inhnity on the left. If —00 < t 2 < 00 , the metric 
( p.5|) covers the upper left half (that is, the left-hand-side exterior and the black hole in¬ 
terior) in the Penrose diagram of this Kruskal manifold. On our spacetime, C 2 covers the 
corresponding regions left of the shell. 


®This last case occurs when e+ = — 1 and e_ = 1 in ( 2.12| ). The spacetime has two bifurcation 
spheres, and the shell passes between them, remaining at all times behind the white-hole and 
black-hole horizons of each infinity |5S]. We are grateful to Eric Poisson for discussions on this 
case. 
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Now, consider our initial data set as a parametrized spacelike liypersurface Sq in this 
spacetime. By the falloff (|2.16|) , Sq is asymptotic at r —cx) to a constant ti hypersurface Si 
in the chart Ci, with r being asymptotic to ri. Similarly, Sq is asymptotic at r —> —cx) to 
a constant t 2 hypersurface S 2 in the chart C 2 , with r being asymptotic to r 2 . Without loss 
of generality, we can take Si and S 2 to be respectively the hypersurfaces fi = 0 and t 2 = 0. 
We assume that Si and S 2 intersect, and that they do so left of the shell, in the black hole 
interior in the left-hand-side Kruskal geometry.^ The value of R at the intersection (where 
R = ri = —r 2 ) is denoted by p. Note that p can be regarded as a piece of gauge-invariant 
information in our initial data set. 

Let So be the hypersurface consisting of Si for ri > p and S 2 for r 2 < —p. So is not 
smooth, but has a corner (a sharp ridge) at ri = —r 2 = p. We choose a positive parameter 7, 
and we take 7 and I so small that (1 -|- 7)p < min(2M_,r — 1 ). We now deform So in the 
regions — (1 + 7)p < 'r 2 < —p and p < ri < (1 -|- 7)/), in a way specihed below, so that the 
deformed hypersurface Sq becomes a smooth, parametrized hypersurface, with a parameter 
r that coincides with ri for r > (l-|-7)p and with r 2 for r < — (l-|-7)p. The canonical data on 
So is by construction gauge equivalent to our original initial data, and it becomes uniquely 
determined after we specify Sq in the transition regions |r| < (1 -|- 'y)p and x — I < r < x. 
As our gauge choice, we take canonical data on So as the representative from the gauge 
equivalence class of our original initial data. 


C. Liouville form and the reduced Hamiltonian theory 


We now hnd the Liouville form Op pulling the Liouville form 6 (|3.3| ) back to the 
transversal surface H. This means that we need to evaluate the right-hand side of (p.3|) 
when the constraints and our gauge condition hold. 

Outside the transition regions |r| < (1 -|- 7 )^ and x — I < r < x, our gauge reads 


R{r) = |r| 
A(r) = 1 , 


(3.6a) 

(3.6b) 


with the gravitational momenta given by (p.2|) . With (p. 6 |) , SR and SA vanish. The only 
contributions to the integral in ( |3.3|) therefore come from the transition regions. We evaluate 
these contributions in appendices ^ and P, specifying the gauge in the transition regions 
and hnally passing to the limit where the parameters I and 7 vanish. From ( [All| ) and ( |B7D , 
we hnd 


= Pp^P + , 


(3.7) 


where 


®This assumption is a further restriction on the initial data. Qualitatively, it tells how “early” or 
“late” the asymptotic ends of Sq may be with respect to each other and the shell trajectory. 
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, /r + p + Vp2 + rn? + J2M,x 
p := v/2iiEr- y2l^ + rln ( - v + 


r + V2M_r 


(3.8a) 

(3.8b) 


with p being a solution to 


M+ — M_ = a/p^ + + —-p 

^ c 


2M. 


(3.9) 


Equation (|3.9| ) has been obtained by eliminating R'_ from ( |B3|) and 

The reduction is thus complete. The functions {p,Pp, t, p) provide a local canonical chart 
on the reduced phase space T, and equations (p.8|) and (|3.9|) determine M+ and M_ as 


(3.10) 


functions in this canonical chart. The Hamiltonian, read off from (|2.17|) , is 

h:=M+ + M_ , 


and the reduced action reads 

As anticipated, T has dimension 4. 


j dt {pfP + pi - ft) 


(3.11) 


D. Dynamics in the reduced theory 


For understanding the dynamical content of the reduced theory, it is useful to introduce 
the new canonical chart (M_,P_,r,p), defined by ( p.8a| ) and 


F_ := 

The new action reads 


— 4Af_ In _ 4^2 M_p . 


S = dt {P_M_ + pir — h) , 


(3.12) 


(3.13) 


where the Hamiltonian h{x,p, M_) is determined by equations ( 3.8b|) -( 3.1(J|) . In this chart, it 
is immediate that both M_ and M+ are constants of motion. It is straightforward to verify 
that the equations of motion for the shell variables are equivalent to (|2.12|) , and thus yield 
the correct dynamics, provided t is identihed as the coordinate H in the spatially flat chart 
( p.4b|) right of the shell. The two solutions of (|T9|) for p correspond to e+ = ±1 in ( p.l2D , 
whereas e_ = 1 always by virtue of the global assumptions made above. We shall provide 
the key steps of this calculation below in section 0 

What remains is the spacetime interpretation of the variable p. Recall that on the initial 
data hypersurface Sq introduced in subsection m. p is the value of R at the sharp ridge 
where the hypersurface fi = 0 in the chart Ci (^.4|), asymptotic to Sq at r —>■ cx), meets the 
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hypersurface ^2 = 0 in the chart C 2 (|3.5|) , asymptotic to Sq at r — —oc. Recall also that 
our Hamiltonian evolves the spacelike hypersurfaces so that at the two inhnities, covered 
respectively by the charts Ci and 6*2, we have dti/dt = 1 and dt 2 /dt = 1. One might 
therefore have thought that as our initial data evolves, p{t) would be the value of R at 
the sharp ridge where the hypersurface ti = t m the chart ( p.4a| ) meets the hypersurface 
0 = t in the chart ( p.5|) . However, this does not hold. The reason is that in the I —>■ 0 
limit, the chart Ci does not reduce to a consistent chart across the shell, not even if one 
were to allow nondifferentiability: the intrinsic metric on the shell history is unambiguous, 
but evaluating this intrinsic metric from the I 0 limit of (3.4a) and from (|3.4b| ) leads 
to mutually inconsistent expressions because the two masses differ. This means that if 
one approaches the shell from the two sides on the “same” constant ti hypersurface, after 
having hrst taken the limit / —0, one arrives at two different two-spheres on the shell 
history. The I 0 limit of one constant ti hypersurface can be interpreted as a continuous 
hypersurface in the spacetime, and this is what we utilized in the gauge choice and the 
evaluation of the Liouville form, but one cannot maintain such an interpretation for a full 
foliation where ti takes values in an open interval. The spacetime interpretation of the 
variable p must therefore be examined more carefully. 

Consider the chart Cf obtained as the I —> 0 limit of the chart Ci left of the shell. 
Denoting the coordinates in Cf by (fi,ri), the metric reads 


ds^ = —dt\ + dri + 


2M_ 

ri 


dti ) + rldVt^, 0 < ri < 


(3.14) 


If r is the proper time along the shell history, we have from ( |3.4b|) and ( |3.14| ) the relation 

2 


= dt\— dr -I- 


2M, 


dti 


= dil — ( dr -I- — dii 


(3.15) 


If we £x the hypersurface = 0 to coincide with the I 0 limit of the initial hypersurface 


fi = 0 for 0 < ri < r, the relation (|3.15|) determines H as a function of H and the shell 
motion, ti = ti{ti). It can now be verihed that p{t) is the value of R at the sharp ridge 
where the hypersurface ti = ti{t) in the chart ( p.l4| ) meets the hypersurface t 2 = t in 
the chart (|3.5|) . The algebra involved in this calculation appears not to be particularly 
instructive, and it will not be reproduced here. 


E. Comments 

As noted above, the details of our reduction relied on certain qualitative assumptions 
abont the shell motion. In particular, we assumed the shell trajectory to intersect the right- 
hand-side exterior region of the Kruskal geometry right of the shell. Our gauge choice, 
involving the ingoing spatially flat coordinates, allows us to follow the shell trajectories into 


16 



















the black hole, but not into the white hole. A time-reversed gauge choice, involving the 
outgoing spatially flat coordinates, would conversely allow us to follow the trajectories into 
the white hole but not into the black hole. 

In the reduced theory (|3.13|) , the value of the canonical coordinate M_ is a constant 
of motion. If we are only interested in the shell motion, we can reduce the theory further 
by dropping the Liouville term P_M_ and regarding M_ as a prescribed positive constant. 
This is arguably the theory of physical interest for an observer who scrutinizes the shell 
motion from one asymptotically flat infinity and regards the “interior” mass as fixed. The 
action then reads 


S 


dt (pf — h) , 


(3.16) 


where h{x,p, M_) is determined by equations (|3.8b|) - (|3.10|) . In the limit m —0, this theory 
reduces to that obtained by Kraus and Wilczek |Q by a less direct Hamiltonian reduction. 


IV. TIME-REPARAMETRIZATION 

In this section we first present a general formalism for reparametrizing time in a Hamil¬ 
tonian system with a two-dimensional phase space. We then apply this formalism to the 
reduced Hamiltonian theory (|3.16|) . 


A. General time-reparametrization formalism for two-dimensional phase space 

Consider a Hamiltonian system with a two-dimensional phase space F := {{q,p)} and a 
time-independent Hamiltonian h{q,p). With respect to a time t, Hamilton’s equations read 


dq 

dt 

dp 

dt 


dh 

dp 

dq 


(4.1a) 

(4.1b) 


We wish to And a Hamiltonian system that generates the equivalent dynamics with respect 
to a new parameter time T, related to t by 


dT = Ndt 


(4.2) 


where V is a prescribed (positive) function of some suitable set of dynamical variables. We 
further wish this time-reparametrization to preserve the value of the Hamiltonian for each 
solution to the equations of motion ([4.1|) . We examine separately two cases: 1) V is a 
function on F, and 2) V is a function of q and the new velocity V := dq/dT. 
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1. N = N{q,p) 


Suppose that N{q,p) is a prescribed function on F. We replace p by a new momentum 
P := P{q,p), where 


^ = /V(,.p) . (4.3) 

We assume that P{q,p) is an invertible function of p for each q, with the inverse p{q,P). 
The new phase space is F := {{q, P)}, and we take the Hamiltonian on F to be 

H{q,P):=hiq,p{q,P)) . (4.4) 

Hamilton’s equations on F with respect to a time T are then easily seen to be equivalent 
to (H). provided t and T are related by (|4.2| ). 


2. N = N{q, V) 


Suppose next that N{q, V) is a prescribed function of q and the new velocity V. 

Recall that equation ([4.1aj) dehnes the velocity v := dq/dt as a function on F. We assume 
that this function can be inverted for the momentum as p = p{q, v). We can then dehne on 
the velocity space the energy function 

h{q,v) := h{q,p{q,v)) . (4.5) 

The dynamics is now encoded in the statement that h{q, v) is constant in t. The value of 
h{q,v) provides one constant of integration, and expressing dq/dt in terms of this constant 
and q yields the general solution in terms of a single quadrature. 

Consider now the time-reparametrization ([4.2|) with N = N{q,V). The velocities v = 
dq/dt and V = dq/dT are related by 


v = N{q,V)V . 

Using (|4.6|) , we can dehne on the new velocity space the energy function 

H{q,V):=hiq,N{q,V)V) . 


(4.6) 


(4.7) 


Provided the relation ( [4.6|) between the velocities is not degenerate, the full 
then encoded in the statement that H{q, V) is constant in T. 

We wish to hnd a Hamiltonian H{q, P) from which P[{q, V) emerges as the 
tion. If L{q, V) is the corresponding Lagrangian, we have 


dynamics is 
energy func- 


and 


H{q,V) 


y dHd, V) 

dV 


Hq,V) 


(4.8) 
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( 4 . 9 ) 

Solving ([4.81) for L{q, V), we find from (|4.9|) that the general solution for P{q, V) is equivalent 
to 


dP(q,V) _ , 8H(q,V) 

8V dV ' ' ’ 

The Hamiltonian H{q,P) is obtained by inverting P{q,V) for V and substituting this in 
Hiq,V). 


3. Comments 


Our time-reparametrization preserves the value of the Hamiltonian on each solution to 
the equations of motion. It does not, however, preserve the value of the action, and it cannot 
in general be thought of as a canonical transformation. 

After N is specihed, the solutions to (|4.3|) and ([4.10|) each contain an arbitrary additive 
function of q. This arbitrariness corresponds to a canonical transformation that redehnes P 
by the addition of (the gradient of) an arbitrary function. 


4- Example: relativistic partiele in (1+1)-dimensional Minkowski space 


As a simple example, we apply this reparametrization formalism to the free relativistic 
particle in (l+l)-dimensional Minkowski space. 

We start from the Hamiltonian 


h = \/p‘^ + 711? , 

which evolves the particle in the Minkowski time t. We then have 

dh p 




dp _|_ ^2 


(4.11) 


(4.12) 


We wish to identify the new time parameter T as the proper time of the particle. From (|4.2|) 
and (|4.12|) we then obtain 


N = Vl -n2 = 


m 




(4.13) 


We can thus use the above formalism with N{q,p) = m (p^ + m?) As a solution to ( [4.3|) , 
we choose P{q,p) = msinh” (p/m). This leads to the familiar point particle proper time 
Hamiltonian 


H(q, P) = m cosh(P/ m) 


(4.14) 
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B. Proper-time Hamiltonian for the self-gravitating shell 


We now apply the time-reparametrization formalism of subsection |IV A| to the Hamilto¬ 
nian theory ( p.l6|) . Our goal is to obtain a Hamiltonian that evolves the shell with respect to 
its proper time. We follow the route of subsubsection [IV A 2| , specifying the reparametriza- 
tion in terms of the new velocity. M_ will be regarded as a prescribed constant throughout. 

We hrst need the Hamiltonian h = M+ + M_ (|3.10 ) as a function of the old velocity f. 
Using the implicit relations ( ^.8b| ) and ( |3.9| ) to evaluate dM^/dp, we hnd that Hamilton’s 
equation f = dh/dp takes the form 


r = 




2M4 


where p is still implicitly given by ( p.9|) . Solving (|4.15|) for p and substituting in 
M^- M_ m 1 - (^ + ^2^+/^) a/2^+A 


m 


m 

Yx 


(4.15) 
yields 

(4.16) 


t -|- 2 A(f_|_ /1 


Equation (|4.16|) determines M+, and hence h, as a function of r and f. 

Let r denote the proper time of the shell. As the parameter time t coincides with the 
spatially flat time ti in the metric (|3.4b|) right of the shell, we have 


dr^ = dv — \dx + 


2M4 


dt 


(4.17) 


This can be solved for dt/dr as 


dt _ A2^+A (dx/dr) + {dx/drf + 1 - 2M+/r 
dr (1 — 2M+/r) 


(4.18) 


where the parameter e+ = ±1 labels the two solutions. Using (|4.18|) to express x in terms 
of dx/dr, we can put equation ([4.16|) in the form 


M+ - M 


m 


m 


= e+y (dx/dr) + 1 — 2M+/r 


(4.19) 


As M_|_ is a constant of motion, the shell motion is completely determined by equation ( |4.19D . 
Comparing ([4.19|) to ( |2.12|) shows that our reduced Hamiltonian theory has correctly repro¬ 
duced the shell motion that arises from Israel’s junction condition formalism, with the pa¬ 
rameter e+ coinciding with the parameter e+ in ( |2.12| ). Equation ( [4.19|) results from squaring 
( P.12| ) once, in a way that eliminates the parameter e_; however, as e_ = 1 by our global 
assumptions, the full information in (|2.12| ) is contained in ( |4.19|) . 

Solving (|4.19|) for M+ yields 


M. - M_ 


m 


+ ^ = ^ (dx/drf + 1 - 2M_/x 


(4.20) 


20 




































































As M_|_ > M_, only the positive sign for the square root in (|4.2CI|) can occur; in terms 
of ( 2.12 ), this sign is equal to e_. From ( 4.20|) , the energy function on the new velocity 
space reads 


H{x, V) = M+(r, V) + M_= m^/V^ + 1 - 2M_/r - — + 2M_ , (4.21) 

2r 


where we have written, in the notation of subsection |1V 7\| , V := dx/dr. As a particular 

(4.22) 


solution to ( 4.10 ) we choose 

P(x, V)=m\n(v+ + 1 - 2M_/r) 


Inverting this for V and substituting in (|4.21|) gives the new Hamiltonian 

,2 

+ M_ 2 - ( 


m 


H{x, P) = mcosh(P/m)-h M_ 

2r 


—) exp(—P/m) 


(4.23) 


C. Minkowski-like Hamiltonian for the self-gravitating shell 


We now consider a time-reparametrization that makes the shell Hamiltonian analogous 
to the Minkowski time point particle Hamiltonian (^4.11| ), which is also the Minkowski time 
Hamiltonian for a free spherical, non-gravitating dust shell in flat space. Starting from the 
shell proper-time Hamiltonian ([4.23|) , we denote the new momentum by p, and we run the 
formalism of subsubsection [IV A 1| backwards with the choice N (r, p) = m (p' + m^) 


As with the point particle example in subsubsection |1VA4| , we solve 


by P(t, p) = 


msinh (p/m). Denoting the counterpart of h in ([4.4]) by h(r, p), we obtain 


h(r, p) = v^p2 m2 - — + M_ 

2x 


^ + rn? 


P 


(4.24) 


V. GEON WITH A SELF-GRAVITATING SHELL 

In this section we adapt the formalism to a shell in a spacetime with the RP^ geon 
topology. 

As mentioned in the Introduction, an asymptotically flat, spherically symmetric space- 
time {M,g) with a single asymptopia can have spatial topology RP^\{a point at inhnity}. 
We refer to an asymptotically flat spacetime with this topology, or to an asymptotically flat 
initial data set in such a spacetime, as an RP^ geon. The covering space of the spacetime 
then has the wormhole topology of the extended Schwarzschild geometry. 

In vacuum, one can obtain a spherically symmetric RP^ geon Einstein spacetime as the 
quotient of Kruskal manifold under a freely and properly discontinuously acting involutive 
isometry [Q. Let (M, g) be Kruskal manifold, and let {t, x, 6, 0) be a chart in which t and x 
are the usual Kruskal time and space coordinates (denoted respectively by v and u in Ref. 
p3|). The isometry in question is then 
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/ : {t, X, 9, 0 ) I— *• (t, —X, 71 — 9, (j) + 7i) 


( 5 . 1 ) 


As I commutes with rotations, the quotient spacetime {M,g) := {M,g)/I is spherically 
symmetric. In {M,g), the constant t hypersurfaces that do not hit a singularity have topol¬ 
ogy X M, with two asymptotically flat inhnities, and they have at a; = 0 a wormhole 
throat at which the radius of the reaches its minimum value. In {M,g), the correspond¬ 
ing constant i hypersurfaces have topology MP^\{a point at inhnity}, and the throat has 
become a “minimum radius” two-surface with topology MP^. Away from the throat history, 
{M,g) is indistinguishable from half (say, x > 0) of {M,g). The Penrose diagram can be 
found in Ref. []^. Note that the throat history in {M,g) is only defined with respect to a 
given foliation, while the throat history in (M, g) has a coordinate invariant meaning as the 
trajectory of the “minimum radius” MP^. The reason for this difference is that / does not 
commute with the Killing time translations on {M,g): these Killing time translations do 
not descend into globally-dehned isometries of (M, g). 

Consider now a spherically symmetric spacetime that has the MP^ geon topology and 
solves Einstein’s equations with a spherical dust shell. Away from the shell, Birkhoff’s theo¬ 
rem still guarantees that the spacetime is locally isometric to Kruskal manifold. We assume 


that the spacetime right of the shell is as in section |T|: this part of the spacetime is part of 
Kruskal geometry, containing the right-hand-side Kruskal inhnity, and the shell trajectory 
intersects the right-hand-side exterior region in this Kruskal geometry. The spacetime left 
of the shell is assumed to be part of the vacuum MP^ geon spacetime described above, and 
to contain the throat history. 

If the shell passes through the throat, it needs to cross itself there. We assume that such 
a crossing does not happen. 

A Cauchy surface in this spacetime has only one infinity, in the part right of the shell, 
whereas the part left of the shell is compact. We can therefore unambiguously regard the 
left-hand side of the shell as the interior and the right-hand side as the exterior. 

It is easy to adapt the Hamiltonian formalism of section to these MP^ boundary 
conditions. We take N, A, R, tta, and tir to be even in r and odd in r, with the 
consequence that M+ = M_ in the falloff (|2.16|) . We assume r > 0, add to the system a 
second shell at r = —r, and hnally take the quotient of the spacetime under the isometry 
(f, r, 6*, 0) H->■ (t, —r, TT — 0, 0 -f vr). The resulting Hamiltonian theory is clearly consistent in 
the same sense as the Kruskal-type theory of section ^ The action can be written as 


S = S^ + S: 




(5.2) 


where Sy, is given by (|2.5|), with the r-integration extending from r = 0 to r = cx), and 


SdT, — 


dt Ml 


(6.3) 


When the equations of motion hold, we recover the above RP^ Einstein spacetimes with a 
dust shell. The throat is located at r = 0. 

The Hamiltonian reduction proceeds in close analogy with that in section B To choose 
the gauge, we introduce the analogue of the chart Ci (|^) , with M_ > 0 now denoting the 
mass in the interior. The range of ri is bounded below by the fi-dependent throat radius. 
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and it is the throat radius that emerges as the parameter p. The transition region near the 
throat is handled as in appendix but because now r > 0 in our action, the contribution to 
the Liouville form is only half of that found in appendix The transition region near the 
shell is handled exactly as in appendix p. The reduced action is given by equations (|3.8|) - 
( p.ll| ), with the exceptions that the right-hand side in the counterpart of (p. 8 a|) contains the 
factor and (p.lO|) is replaced by 


h := M. 


(5.4) 


From section |T| it is clear that the reduced theory reproduces the correct equations of 
motion. In the classical solutions, p{t) is the value of R at the throat in a foliation dehned 
as with the chart (|3.14|) . 

A canonical transformation that replaces the pair {p,Pp) by leads to the ac¬ 

tion (|3.18| ), with ( |3.1U|) replaced by (|5.4| ). Dropping the term P_M_ gives a theory in 
which the interior mass M_ is regarded as a prescribed positive constant. The time- 
reparametrizations of section 0 clearly carry through without change: the counterparts 
of the Hamiltonians (|4.23|) and (|4.24|) differ only in that the (constant) additive term 2M_ 
is replaced by M_. 


VI. SELF-GRAVITATING SHELL WITH SPATIAL TOPOLOGY 


In this section we consider the spatial topology 

We start directly from the action principle. In the bulk action (|2.5|) , we take 0 < r < oo, 
with the falloff ( p.l 6 | ) at r 
introduce the falloff 


—> cx). The total action is given by (p.2|) and (|0|). 

At r —>• 0, we 

A{t,r) = Ao + O{r‘^) , 

(6.1a) 

R{t,r) = Rir + 0{r^) , 

(6.1b) 

7rA(t,r) = TTAar^ + O(r^) , 

(6.1c) 

7iR{t,r) = TiR^r + 0{r^) , 

(6.Id) 

N{t,r) = No +, 

(6.1e) 

W(t,r) = V[r+ 0(r3) , 

(6. It) 

2 , tir^, No > 0, and iV[ are functions of t only. It 

is straightfor- 


where Aq > 0, i?i > 0, ttaj 
ward to verify that this falloff is consistent with the constraints and preserved by the time 
evolution, and no additional boundary terms in the action are needed at r = 0. From (p.2|) 
we see that in the classical solutions, the mass left of the shell must vanish, and r = 0 is 
just the coordinate singularity at the center of hyperspherical coordinates in flat space. The 
classical solutions therefore describe a self-gravitating shell with a flat interior. The spatial 
topology is M^. 

The reduction proceeds as above, using the analogue of the chart Ci (|3.4|) with M_ = 0 
and ri > 0. In the region ri < r — the initial data hypersurface Sq extends smoothly 


to ri = 0, and there is no counterpart of the parameter p of the Kruskal and MP'^-geon 
topologies. The only contribution to the Liouville form comes from the shell transition 
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region, which is handled exactly as in appendix ^ but with M_ = 0. The reduced action is 
given by (|3.16|) and (|^) , where M+ is obtained from (|3.8b|) and (|3^) with M_ = 0. It is 
again clear from section 0 that this reduced theory reproduces the correct dynamics. As 
expected, the reduced phase space is two-dimensional. 

The time-reparametrizations of section carry through without change. The coun¬ 
terparts of the Hamiltonians (|4.23|) and ( [f.24|) are obtained from these formulas by simply 
setting M_ = 0. In particular, (|4.24|) reduces to the Hamiltonian used in Refs. 


VII. REMARKS ON QUANTIZATION 

In this section we discuss the prospects for quantizing the reduced theories. We first 
review the pure vacuum case, and then turn to the coupled system. 


A. Mass spectrum of spherically symmetric vacuum wormholes and MP^ geons 


In sections EHYII we considered the dynamics of a shell coupled to spacetime geometry. 
However, the methods immediately adapt to spherically symmetric vacuum gravity by simply 
omitting the shell. 

With the Kruskal topology, and no asymptotic masses fixed, the constraints imply M+ = 
M_ := M. With the gauge choice of section |T|, without the shell, the reduced action reads 


S = dt {ppp - h) , 


(7.1) 


where h = 2M, and M is obtained from (|3.8a|) with M_ = M. Geometrically, p{t) is the 
value of R at the sharp ridge in the foliation described in subsection Em without the shell; 
this ridge evolves in the black hole interior along a history of constant Killing time. At 
t —>■ —oo, we have p{t) 2M as the ridge approaches the bifurcation two-sphere, but in 
the future the gauge breaks down at a finite value of t as p{t) —> 0. With the MP^-geon 
topology, the only differences are that h = M, and the right-hand side of ( |3.8a| ) contains 
an additional factor i. The reduced phase space is two-dimensional in each case. These 
results agree with those obtained by Kuchaf’s reduction method [^,0 under a falloff that is 
qualitatively similar but makes the constant t hypersurfaces asymptotic to hypersurfaces of 
constant Minkowski time. 

If one chooses to fix the mass at one infinity with the Kruskal topology, the reduced 
theory has no degrees of freedom. The same holds if one chooses to fix the mass at the 
infinity or at the throat with the MP^ geon topology. With the topology, the reduced 
theory is always void. 

Quantizing the reduced theories with a zero-dimensional reduced phase space is of 
course trivial: the mass M is a prescribed c-number. Quantizing the theories with a two- 
dimensional reduced phases space offers, however, several options. 

One option is to perform first a canonical transformation to the pair (M, Pm) as 
in sections |l| and 0. One can then take quantum states to be described by func¬ 
tions ^(M) of the positive-valued configuration variable M, adopt the inner product 
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^$i|^ 2 ) = /g°° (iM ^i(M)^ 2 (Af) (or a similar inner product with some M-dependent weight 
factor), and promote M into the quantum operator M that acts in the Schrodinger picture 
as IHi 




(7.2) 


The spectrum of M, and thus also that of the Hamiltonian operator h, is continuous and 
consists of the positive real axis. 

Another option is to take quantum states to be described by functions V’(p) of the 
positive-valued “throat radius” p, adopt the inner product {'ipi\'ip 2 ) = f^{p)dp 4 ’i{p)'ip 2 ip) 
where p(p) is some weight factor, and try to promote the function M{p,pp) into an operator 
on this Hilbert space. As our M{p,pp) is known only implicitly, we have not tried to pursue 
this quantization, but there seems no obvious reason to expect that the spectral properties 
of the resulting Hamiltonian operator would agree with those of the operator M in (|7.2|) . 

Indeed, quantization of spherically symmetric vacuum gravity was discussed in Ref. 
in terms of a related “wormhole throat” phase space {a,pa), on which the Schwarzschild 
mass is given by 


M(a,pa) = I (^ + a 


(7.3) 


The conhguration variable a has an interpretation as the radius of the wormhole throat, 
much as our p, but with a time parameter that is now identihed with the proper time of the 
throat history.]^ If the Hilbert space is chosen as above with the conhguration variable p, 
with reasonable choices for the weight factor p, the function M{a,Pa) can be promoted into 
a self-adjoint operator whose spectrum is bounded below and purely discrete [|^. 

We regard as artihcial the continuous mass spectrum arising from the quantiza¬ 
tion (Q). because one can similarly obtain a continuous spectrum for any dynamical sys¬ 
tem whose Hamiltonian is not explicitly time-dependent. For any function H with non¬ 
vanishing gradient on the phase space, one can hnd a local canonical chart of the form 
{H, ^ 2 , • ■ ■) Qn,PH,P 2 , ■ • ■ ,Pn), in whlch H is one of the canonical coordinates. If the range 
of H in this chart is M+, one can adopt a Schrodinger representation with Hilbert space 
L2(R+)<»H, with H a Hilbert space for the remaining g’s. The Hamiltonian operator H can 
then be taken to act as a multiplication operator, 

H'il!{H,q 2 ,-■ ■ ,qn) = H'4!{H,q2,-■ ■ ,qn) , (7.4) 


and its spectrum is M+. 


^^The Hamiltonian M{a,pa) (|7.3D describing the proper-time evolution of the throat was pre¬ 
viously considered by Friedman, Redmount and Winters-Hilt [^J60t| without a derivation by re¬ 


duction from spherically symmetric vacuum gravity. In Ref. |]^|, this Hamiltonian was derived 
from Kuchaf’s reduced Hamiltonian theory Q by a canonical transformation. A similar derivation 
could clearly be given from the canonical pair (M, Pm) of the present paper, despite the technical 
differences in our falloff and that of Ref. [Hi . 
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Ambiguities in canonical quantization are, of course, well recognized |]6l|-|63 


One specific 

issue not addressed above is in the global properties of the canonical transformations. For 
example, the canonical transformation that takes the phase space (a, Pa) to Kuchaf’s reduced 
phase space 0 is not onto: the classical dynamics in Kuchaf’s reduced phase space is 
complete, but the classical dynamics in the phase space (a, pa) is not One’s attitude to 


such classical incompleteness in view of quantization may depend on what one sees as the 
role of singularities in quantum gravity P!^,|5^-[70 . 


B. Quantization of shell coupled to geometry 


We now turn to the coupled system. We restrict consideration to the proper-time Hamil¬ 
tonian ([4.23| ). 

When M_ = 0, the shell encloses a flat interior with trivial topology, and the Hamiltonian 
([4.23|) takes the form corresponding to a relativistic particle in a Coulomb potential. 


= mcosh{P/m) — 


m 


(7.5) 


discussed by Hajicek One can adopt a Schrodinger representation corresponding to 

conhguration-space variable r G M+ and Hilbert space 


H ;= L2(M+,r"dr) 

where a is a parameter. With the factor ordering 


N 


[cosh(P/m)] = lim r 

N —vrvn f ^ 


■ 1 )' 


N^oo (2n)\ 


,^n^a/2 


n=0 


where A = di, H becomes a self-adjoint operator H with domain 


D{H) = {/ I /'"">(0) = 0, /!“> € L 2 , all n} 


(7.6) 


(7.7) 


(7.8) 


(Hajicek takes a = —1, but notes the unitary equivalence of for a different choice 

of a.) For m < 1.9, H is bounded below, and its spectrum, like that of the nonrelativistic 
Coulomb problem, has discrete and continuous parts. 

When M_ > 0, one expects that the Hamiltonian (|4.23|) can be made into a self-adjoint 
operator in an analogous manner, and one expects the spectrum then to be bounded below 
and partly discrete for small values of m. However, there appears to be no reason to expect 
that the term proportional to M_ would allow the spectrum to have a lower bound for large 
values of m. Oharu and Winters-Hilt 0 are currently examining a self-adjoint extension 
of H on L 2 (M+, dr), with factor ordering corresponding to the choice (|7.7|) with a = 0: 

2 

H = m[cosh(P/m)] — — -M_mr“^'’'^[exp(—P/m)]r“^/^ -|- 2M_ . (7.9) 

2r 

Finally, recall that our time-reparametrization derivation of the proper-time Hamiltonian 
( |4.23| ) assumed M_ to be a prescribed, time-independent constant. With spatial topol¬ 
ogy this assumption is automatically satished. With the Kruskal and MP^-geon topologies. 
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on the other hand, one could ask whether it is still possible to carry out an analogous time- 
reparametrization when M_ is a dynamical variable and the phase space is four-dimensional. 
If the answer is affirmative, one could presumably raise anew the issues regarding the spec¬ 
trum of M_ that were addressed in the context of the vacuum theory in subsection |VII A . 
If, after the reparametrization, the dynamics of M_ still decouples from the dynamics of the 
shell as in sections 0 and 0, one could effectively separate variables by hrst considering the 
eigenvalue equation for M_, 


M-ip = M^xjj 


(7.10) 


For each eigenspace of M_, the shell Hamiltonian would then have the form (|4.23|) with 
a c-number M_, and the character of the total spectrum would depend on the spectrum 
of M_. 


VIII. SUMMARY AND DISCUSSION 

In this paper we have considered the Hamiltonian dynamics of spherically symmetric 
spacetimes that contain an idealized, inhnitesimally thin massive dust shell. We considered 
the Kruskal-like spatial topology x M, the MP^-geon spatial topology MP^\{a point at 
inhnity}, and the Euclidean spatial topology The variational equations that arose from 
the unreduced Hamiltonian action were not strictly consistent in a distributional sense, but 
we were able to localize the ambiguity into the single equation that arises by varying the 
action with respect to the shell position. When the ambiguous contribution to this equation 
was interpreted as the average of its values on the two sides of the shell, we correctly 
reproduced the content of Israel’s junction condition formalism. 

We performed a Hamiltonian reduction by adopting a gauge with piecewise flat spatial 
sections, and passing to the limit in which the interpolating transition regions became van¬ 
ishingly narrow. The constraints could then be explicitly solved. For the Kruskal and MP^ 
topologies the reduced phase space was four-dimensional, with one canonical pair closely 
associated with the shell motion and the other pair with the dynamics of the geometry. 
In the limit where the shell is not present, this correctly reproduced previous results for 
spherically symmetric vacuum geometries. Retaining the shell but prescribing by hand one 
asymptotic mass for the Kruskal topology, and the interior mass for the MP^ topology, we 
recovered theories whose reduced phase space was two-dimensional, with just the canoni¬ 
cal pair associated with the shell motion surviving. For the topology, the interior mass 
necessarily vanishes, and we only obtained a two-dimensional phase space, with the single 
canonical pair describing the shell motion. 

For each of the three spatial topologies, we time-reparametrized the dynamics in the 
two-dimensional phase space that describes the shell motion with hxed interior mass. With 
one choice for the reparametrization, we recovered a previously known Hamiltonian that 
evolves the shell with respect to its proper time. With another choice, we recovered a 
Hamiltonian analogous to the square-root Hamiltonian of a spherical test shell in Minkowski 
space. Finally, we briefly discussed the spectra that would be expected to emerge in different 
approaches of canonically quantizing the theories. 
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Our results provide a robust description of the reduced Hamiltonian dynamics of a spher¬ 
ically symmetric dust shell coupled to gravity, in the region of the reduced phase space that 
is covered by our piecewise spatially flat gauge. While this gauge is not global, one can ar¬ 
gue that this gauge and its time-inverted counterpart cover the region of the reduced phase 
space that is of interest to an observer who scrutinizes the shell motion from one asymp¬ 
totically flat inhnity. What remains open, however, is the global structure of the reduced 
phase space. One would also like to describe the reduced phase space in a way that is more 
geometrical and less tied to a particular gauge. One possible avenue for this, currently under 
investigation by Hajicek and Kijowski BUI, might be to generalize to the massive shell the 
canonical transformations that Kuchaf introduced to simplify the vacuum theory |p. Work 
on the analogous problem with a null-dust shell is in progress [Q . 


More ambitiously, one would like to consider systems with matter that is more interesting 
than a dust shell. The canonical formulation of Einstein gravity coupled to a continuous 
distribution of massive or null dust has been discussed respectively in Refs. and |]7n| . 
For the canonical formulation in the presence of other types of fluids, see for example Ref. 
0 and the references therein. A discussion of the difficulties involved with spherically 
symmetric gravity coupled to a scalar held is given in Refs. |78|J79[| . A discussion in the 
context of a dilatonic black hole can be found in Ref. [BD . 
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APPENDIX A: RIDGE TRANSITION REGION 

In this appendix we specify the gauge in the ridge transition region |r| < (1 -|- yjp, and 
evaluate the contribution from this region to the integral on the right-hand side of (|3.3|) in 
the limit 7 — 0 . 


1. Gauge choice 


To specify the gauge in the region |r| < (1 -|- yjp, we consider the classical spacetime of 
subsection Eng. and the spacelike hypersurface Sq in this spacetime. The part |r| < (Ij-yjp 
of Eq li^s in the black-hole region of the Kruskal spacetime left of the shell. 

Let h : M —> R be a smooth function such that 


h{x) 


0 , a: < 0, 

X - I , 2 : > 1 , 


and cPh/dx^ > 0 for 0 < x < 1. We write h^^\x) := d'^h{x)/dx"-. 
For |r| < (1 -|- 7 )p, we seek a gauge in the form 


(Al) 
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R{r) = 


\r\ — p 
IP 


7M^)+l + h 


+ ^0 ) 
P , 


(A 2 a) 

(A 2 b) 


where Aq is a positive parameter. In the subregion |r| < p, the radius of the two-sphere 
is constant on Sq, i? = (1 -|- ^ 7 )p, and the proper distance on Sq is Aodr. The subregions 
P < kl < (l + 7 )p interpolate smoothly between this constant radius gauge and the spatially 
flat gauge (p. 6 |). Note that A(r) > 0, and (1 -|- \'))p < .R < (1 + 7 )p. 

Recall from subsection p!IIB| that (1 -|- 7 )p < 2M_. Equations ( p.2|) thus yield a real¬ 
valued solution for tta and ttr for all of |r| < (1 -|- 7 )p. The gauge ([A^ ) therefore specihes 
a spacelike hypersurface in an interior Kruskal geometry with mass M_, with the ends at 
R = (l-|- 7 )p. What remains is to choose the parameter Aq in (|A2a|) so that this hypersurface 
precisely hts between the points |r| = (1 -|- 7 )p. 

In the curvature coordinates (T, R) in the black hole interior, the metric reads 


ds^ = — 


2M_ 


-1 

- 1 1 dR^ + 


2M_ 


-l\dT^ + R^dQ^ 


(A3) 


where 0 < i? < 2M_, R decreases to the future, and we take T to increase to the right. The 
transformation from (|A3|) to the chart Ci of subsection [III B| reads 


T = ti — 2\/2M_ri — 2M_ In ( ^ ) -|- constant , (A4a) 


\ y/2M_ -\- yAq 


R = ri 


(A4b) 


As our (prospective) deformation ( |A^) of Sq to Sq is symmetric around r = 0 , the value of 
T at r = 0 on Eq is the same as the value of T at the unsmoothed ridge on Sq- On Eq, we 
thus have 


Tr=(i+ 7 )p ~ Tr=Q — 2 '\/ 2 M_p — 2 '\/ 2 M_(l -|- 7 )p 

(V 2 M — -s/p) (y/2MZ + \/(l + 7 )p 


+ 2 M_ In 


On the other hand, from Eq. (80) of Ref. [^] we have 

Atta 


{V2M_ + y/p) f ^/2MZ — 1/(1 + 7)p 


(AS) 


T' = 


2M_-R 


(A 6 ) 


Integrating ( |A 6 |) from r = 0 to r = (l-|- 7 )p, with tta given by (|3.2a| ), and equating the result 
to ([Asl), gives a relation that implicitly determines Aq in terms of M_, p, and 7 . By the 
symmetry of Eq around r = 0 , the relation obtained by similarly comparing to 

rr=o, using the chart 02 , contains exactly the same information. This completes the gauge 
choice. 

We shall below be interested in the limit of small 7 . In this limit, the relation determining 
Aq admits a power series expansion in 7 . The result is 
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7 


An — 


+ 0 ( 7 ^) , 


(A 7 ) 


2^l-p/{2M) 

where O stands for a 7 -dependent quantity that is bounded by a constant times its argument. 


2. Liouville form 


r| < 


We now evaluate the contribution to the integral in the Liouville form (|3.3|) from 
(1 -|- 7 )p, in the limit 7 — 0 . 

As we have noted, the gauge (|A 2 |) for |r| < (1 -|- 7 )/) joins smoothly to the spatially 
flat gauge outside this interval, and the expressions given in are in fact valid for all 
of —cx) < r < X — 1. The differentials (5A(r) and dR{r) therefore contain no delta-functions 
in r at |r| = (1 -|- 7 )^, and it is sufficient to consider the contributions from |r| < p and 
p < |r| < (1 - 1 - 7 )p. 

For |r| < p, we have h = = 0. Equations 


SR = 0 ( 1 ), and equations 
is therefore 0(7). 

Suppose then p < r < (1 -|- 7)p. We now obtain 


and ([A 7 D yield SA = 0(7) and 
yield tta = 0 ( 1 ) and tir = 0(7). The contribution to 


SA = - 6p + 0(l) , 

IP 

SR= (l-hA)) 5 p + 0 ( 7 ) , 
tta = P^/{R'/AY - 1 + 2M_/R + 0(7) , 

, Pin 

^/{R'/Af - 1 + 2M_/R 

where the argument of h and its derivatives is always (r — p)Hp. Note that the hrst term 
in i\KEdi ) is 0(7-1). 

For J drn\SA, changing the integration variable from r to a; := (r — p)/7p gives 

/‘{l+7)P /•! 

/ dr7iAdA = —pSp / dxh^‘^\x)\/{R'/AY — 1 -|- 2M_/R -1-0(7) 

J p Jo 

= —^j 2 M_p 5 p f dxh^‘^\x) -I-0(1) 

Jo 

= -^2M_p5p + 0 ( 1 ) , (A9) 


(A8a) 

(A8b) 

(A8c) 

(A8d) 


where o(l) stands for a 7 -dependent quantity that goes to zero as 7 —0. We have used 
the fact that '^/((R'TAO^'^^T+'TMLT^R —■>/ 2M _/p pointwise in a: as 7 —0, and taken the 
limit under the integral by dominated convergence. 

For f dr tirSR^ the assumption > 0 allows us to change the integration variable from 
r to M := R'/A. We obtain 


»(l+7)p 


dr tirSR 


pSp 


»(l+7)p 


dr (i?VA)'(1 - hA)) 

^/iWjAj^^^TT2Ml/R 


+ 0(7) 
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= p^p 
= p^p 



du (l — 


lo V(i?VA)2 - 1 + 2M_/R 
du 


+ 0(7) 


lo a/m^ - 1 + 2M_/p 


+ 0(1) 


Vv2wr-v^r 


(AlO) 


We have used the facts that -^Z(-R'/A)^ — 1 + 2M_/R a/m^ — 1 + 2M_/p and = 
mAq [1 — m( 1 — Ao)]~^ —> 0 pointwise in u as 7 —>• 0 , and taken the limit under the integral 
by dominated convergence. 

Adding the identical contributions from the region —(l + 7 )p < r < —p, we find that the 
total contribution to the Liouville form (|3.3| ) from the ridge transition region |r| < (1 + 7 )p 
is 


»(l+7)p 


-{1+7)P 


dr ( 7 rA( 5 A + TipidR) = 


In / V2M +Vp \ _2^2M_p 

^ \V^-Vp) ^ 


(5p + o(l) . (All) 


APPENDIX B: SHELL TRANSITION REGION 


In this appendix we specify the gauge in the shell transition region x — I < r < x, and 
evaluate the contribution from this region to the integral on the right-hand side of (|3.3|) in 
the limit I —>■ 0. 


1. Gauge choice 


To specify the gauge in the region x — l<r<x,we again consider the classical spacetime 
of subsection Eng. and the spacelike hypersurface Sq in this spacetime. The part x — I < 
r < r of So lies in the Kruskal spacetime left of the shell. 

Let / : M —>• R be defined by 


/(^) 


0 , X ^ {0, 1). 


(Bl) 


We write •= d^f{x)/dx'^. f is continuous everywhere, and smooth except at x = 0, 

with f^^\x) —>• 1 as X —> 0+ and f^^\x) —> 0 as x —>• 0_. Note that f^‘^\x) —> 0 as x —> O-t. 
For (l-|- 7 )p<r < 00 , we choose the gauge 


A = 1 , 

^ l\/p’^~Trn^ „ f X - r\ 

^ = ^ - ^ 


(B2a) 

(B2b) 


Outside the shell transition region x — I < r < x, this clearly agrees with the spatially flat 
gauge (p.6|) . To show that the gauge is admissible, we hrst note that for r — / < r < r. 
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the constraints are solved by the gravitational momenta given by (|3.2|) . At the shell, the 
Hamiltonian constraint (p.l3a|) is identically satished. The momentum constraint ( p.l3b ) at 
the shell reads, using ( |3.2a)) . 


p = t J {R'_f - 1 + 2M_/r - V2M+r , 


where (B2b) gives 


R'_ = l + 


a/P^ + 


(B3) 


(B4) 


The constraints can therefore be solved both at the shell and away from the shell, and the 
gauge is thus admissible. The gauge is smooth everywhere except at the shell, and at the 
shell it is consistent with the regularity assumptions of section 


2. Liouville form 


We now evaluate the contribution to the integral in the Liouville form ( p.3|) from x — l< 
r < r, in the limit 1^0. 

As the gauge (|B2|) is smooth for (1 + 7 )p < r < cx) except at the shell, the differentials 
<5A(r) and 6R{r) do not contain delta-functions in r except possibly at r = r. Equation 
shows that <5A(r) = 0 everywhere. It is therefore sufficient to consider separately 


ttrSR for r — / < r < r, and the delta-function contribution to TrjidR at r = x. 
For X — I < r < X, (|B2b[ ) gives 


R' = 1 


a/P^ -I- 




R" = 


\/p2 -|- 
lx 




dR = {l - R')St + 0 (l) , 

where the argument of / and its derivatives is (r — r)/I. From (|3.2b| ) we have 

xR" 


ttr = 


0 ( 1 ) , 


(B5a) 

(B5b) 

(B5c) 


(B6) 


^/W^^YT2MIJx 

where we have used the observations i?" = 0{l~^) and R = x + 0{l). Note that the hrst 
term in ( |B6|) is We thus obtain, changing the integration variable from r to v := R', 

drR''(l-R') 


.X 

/ dr7rjiSR = xdx 

Ix-i 


= r5r 


x-i ^/R'^ - 1 + 2M_/x 
dv (1 — v) 


+ 0(0 
+ 0(0 


1 ^/^R^^T+^MZJx 

\/2MIx — \j2M+x — p 

r + p + + A/2M+r 


+rln 


r+ ^j2Mlx 


dx + 0(0 , 


(B7) 
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where we have used (|B4|) for R'_. 


What remains is the delta-function in 6R at r = r. From ( B2b|) we have 




6{r — x) + (non-distributional function of r) . (B8) 


X 


From (|3.2b|), we have 



V2M+/t , 


(B9a) 


(p + v/2iW+r)^ - M_x 


(B9b) 


r(p + yiM+r) 


where we have used ( |B4p and the fact that R'!_ = 0. As tir is not continuous at r = r, the 
product 'KrSR is not dehned as a distribution, and the contribution to the Liouville form 
is ambiguous. However, as are both 0(1), and the delta-function in 5R (P8|) is 0(/), 
we argue that the ambiguous contribution can be taken to vanish in the limit I —>■ 0. It is 
seen in the main text that this leads to a reduced Hamiltonian system that reproduces the 
correct dynamics. 

The ambiguity in tt^SR appears to have the same origin as the ambiguity of the equation 
of motion (|2.81]) in the unreduced formalism: both involve varying the action with respect 
to X. Note that if the function / had been chosen so that 0 as a: —> 0+, R'!_ and vr^ 

would be nonvanishing and proportional to and the above argument for the vanishing 
of the ambiguity in the limit I —> 0 would not apply. 
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